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Abstract 



A new method, analogous to super symmetry transformation in quantum 

mechanics, is developed for a particle in the lowest Landau level moving in 

an arbitrary potential. The method is applied to two-dimensional potentials 

formed by Dirac delta scattering centers. In the periodic case, the problem 

is solved exactly for rational values of the magnetic flux (in units of flux 

quantum) per unit cell. The spectrum is found to be self-similar, resembling 

the Hofstadter butterfly Q. 
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In recent years, the energy spectrum of two-dimensional (2D) electron systems have 
attracted great interest, because of the relevance of the problem to the magneto-transport 
properties and in particular to the quantum Hall effect [|J . It is believed that the physics 
of the integer quantum Hall effect is governed by the interaction of electrons with a disordered 
potential which leads to localization of the eigenfunctions. Observed conductance steps can 
be explained by a sequence of localization-delocalization transitions. Quantization of the 
Hall conductance due to periodic potentials has also been studied fH|J. In this case, the 
presence of steps is explained by the gaps in the energy spectrum. 

There are two opposite approaches to the problem of an electron moving in a periodic 
potential: The tight-binding and the nearly-free-electron methods. In the first approach, the 
magnetic field is introduced via Peierls substitution where the matrix elements are multiplied 
by exp[|^ j A • dl] 0. On the other hand, in the case of nearly- free-electrons, Landau level 
structure is essential and the lattice potential is introduced via intra- and inter-Landau 
level scattering matrix elements. The duality between the position and the momentum 
in quantum mechanics leads to similarities between the two methods. In the presence of 
magnetic field, the secular equations for the two limits, the tight-binding and the nearly- 
free-electron approach, with certain approximations, are identical ||. The characteristic 
feature of the problem is that the secular determinant for the limit of infinite crystals can be 
reduced to a finite determinant, when the magnetic flux per unit cell is a rational number in 
units of the flux quantum. Tight-binding case has been studied by Azbel' and Hofstadter 
]l||, who showed that the system has a complicated self-similar spectrum. 

In this work, a new approach is developed for a particle in the lowest Landau level 
moving in an arbitrary potential. The method is used to obtain the energy spectrum in 
the presence of Dirac delta potentials. The difficulty with Dirac delta potential is that, 
even if the inter-Landau level couplings are small and can be neglected, there is a strong 
intra-Landau level mixing. The formalism developed in this study, leads to an eigenvalue 
problem where the coupling between distant sites become negligibly small. Moreover, since 
the problem is formulated in real space, the distribution of Dirac delta potentials can be 

2 



arbitrary. In spite of the sharpness of the potential, the assumption that particle is confined 
into a single Landau level is consistent with the result that shifts in energy eigenvalues are 
very small in comparison to inter-Landau level separation. The method can be used for any 
potential, as far as the particle is confined into the lowest Landau level. After presenting the 
method, the case of Dirac delta scattering centers will be solved for periodic distribution. 

The Hamiltonian for a particle of mass m and charge q, moving in 2D in the presence 
of magnetic field B = V x A perpendicular to the plane and potential V, is given by 
H = H + V where 

*. = 5j;<p-?A)» . (i) 



Using the symmetric gauge A = |Bxr and complex coordinates z = X+iY = y qB /2hc(x+ 
iy) where r = (x,y), the unperturbed Hamiltonian can be written as H = a) a where 
a* = —d/dz + z*/2. Since [a, a*] = 1, the energy eigenvalues are given by E n = huj(n + 1/2) 
where u> = qB/2mc (q is assumed to be positive) and n— 0, 1, 2, .... The ground state wave 
functions are of the form f(z)e~\ z \ 2 / 2 . Here f(z) is any analytic function of z. A basis for 
the infinitely degenerate ground state is formed by 

<f> m (z,z*) = ^=e-^ 2 (2) 
y-nml 

with m = 0,1,2,.... When the magnetic field is very high the particle is confined into 
the lowest Landau level. This is a good approximation as long as the potential is small in 
comparison to Landau level splitting Hcu. The problem is reduced to the diagonalization of 
the matrix V whose entries are given by 

<m\V\m'>= jj 4>* m {z,z*)V{z,z*)4> m {z,z*)d?z , (3) 

where d 2 z = dxdy. Let A and B be defined by 

A m , zz * = <l>* m (z,z*)yJV(z,z*) (4) 



B zz * tm = y/V{z,Z*)^) m (z,Z*) . 
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Therefore V = AB. For positive V(z,z*), A = B\ Now consider the operator V = BA. 
Any eigenvalue E associated with the eigenket \if)> of the operator V = BA is also an 
eigenvalue of the operator V = AB, except when A\ip>— 0. The corresponding ket is given 
by \ip>— E~ l l 2 'A\ip>. Similarly, any eigenvalue E associated with the eigenket \ip> of the 
operator V = AB is also an eigenvalue of the operator V = BA, except when B\tp>— and 
the corresponding ket is given by \i)>= E~ l / 2 B\i)>. Therefore it is enough to diagonalize 
V , instead of V , to find the energy eigenvalues, except for the ground state. This method is 
analogous to supersymmetric quantum mechanics pioneered by Witten ||. Substitution of 
A and B yields 



<zz*\V\z'z'*>= ^Y^^-WW^Y^l . (5) 

When V(z, z*) is negative V = BA is not Hermitian but the eigenvalues are still real. The 
method can be applied to any problem in the lowest Landau level but its advantages will be 
demonstrated for the potential 

V(r) = V £<%-!•,) . (6) 

i 

In this case, z and z' become discrete variables and the matrix to be diagonalized is 

<i\V\3>= ^ e **H*l a /2-W a /2 . (7) 

7T 

It is interesting to note that the matrix elements are very similar to those of the tight-binding 
problem. Phase is nothing but the usual Peierls factor. <z|V r |j> and <m\V\m'> have the 
same spectrum except for the ground state. Let Ni be the number of impurities and N be 
the number of states both of which go to infinity so that Ni/N — * p. Hence V and V have iVj 
and N eigenvalues, respectively. The bigger one has |iVj — N\ zero eigenvalues. If Ni < N, 
N — Ni eigenvalues of V will be zero. V filters the nonzero eigenvalues of V. However, if 
V itself gives zero eigenvalues, the situation must be examined carefully to distinguish the 
zero eigenvalues from the extra zeros. Once V is diagonalized, the eigenfunctions of V can 
be constructed by acting A on the eigenvectors of V. The impurities can be distributed in 
an arbitrary way. 



The problem of an electron in a magnetic field interacting with point impurities has been 



discussed extensively in the literature [|IO|-|T5|j . It has been shown that the zeros of the wave 



function can be adjusted to coincide with the locations of the scatterers if the concentration 

of the scatterers is low enough. This corresponds to Ni < N case discussed above. The 

presence of N — iVj energy eigenvalues unaffected by the impurities is a direct consequence 

of the dimensionality of the matrix to be diagonalized. Furthermore, the method gives the 

full energy spectrum for any distribution of impurities. It must be noted that Dirac delta 

potentials have some peculiarities in 2D. They do not cause scattering, irrespective of their 

strength. However, in the presence of a magnetic field new energy and length scales are 

introduced, and they lead to nontrivial spectra. 

For regularly distributed impurities the problem can be handled by introducing periodic 

solutions. For simplicity the square lattice case will be investigated. Generalization to an 

arbitrary lattice is straightforward. Now, the impurity positions are given by r mn = d(m, n) 

where m and n are integers and d is the lattice constant. The matrix elements become 

, ~ , / / V r Am - m') 2 + (n - n') 2 . . 

<mn\V\mn>= — exp[— 7r — ( \- i(rnn — nm ))\ (8) 
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where is the magnetic flux passing through a unit cell (or flux per impurity) and 0o is the 
flux quantum (note that the ratio of the area of the system to the number of states is n). 
When 0/0o is a rational number the system becomes periodic. Let 0/0o = a = p/q where p 
and q are relatively prime integers then <mn\V\m'n'> is left invariant under 2q translations 
along lattice axes. If p is even, the matrix elements become periodic with period q instead 
of 2q. In any case, as long as a is rational, the eigenvectors \nk> have the property that 

<Rr|nk>=<0r|nk> e* R (9) 

where R fixes the 2q x 2q (or q x q) unit cell and r denotes the position of the impurity in 
the cell. Here n is an integer (band index) and k is a two dimensional vector. Hence, the 
problem is reduced to determination of the coefficients <0r|nk> which satisfy 

Kr'(k) <0r>k>= E nk <0r|nk> (10) 



where 



V TT ,{k) = <Rt\V\Ot'> e 



(11) 



R 



Given the flux ratio a, the energy eigenvalues _E„k can be found by diagonalizing 4q 2 x Aq 2 
matrix V^v(k). However, this does not imply that there will be 4q 2 energy bands. If the 
problem is formulated in the Landau gauge, the additional symmetry can be seen easily. It 
turns out that unit cell has dimensions 1 x q and hence there will be q bands. The symmetric 
gauge leads to a simple formulation of the problem but in the mean time some symmetries 
are obscured. 

In Fig. m energy bands for a = 4/3 and a = 10/7 are plotted. There are 3 and 7 
energy bands, respectively. Using the relation \ip>= E~ 1 / 2 A\ip>, it is possible to find the 
corresponding eigenfunction ip(r) of the original problem. Since ijj(r) belongs to the lowest 
Landau level it must be in the form of f{z)e^ z ^ I 2 where f(z) is any analytic function. 
Analytic functions are determined by their zeros. Considering the translational invariance 
of the system one expects to have regularly distributed zeros. Figure [2] shows the contour 
plots of I^K 1 ")! l° r a = 2- This value of the flux ratio is very special in that all of the matrix 
elements <mn\V\m'n'> become real and the problem reduces to a tight-binding system. For 
this particular case, energy eigenvalues and eigenfunctions can be found explicitly in 
terms of the theta function 



which vanishes at z = (I + 1/2) + (m + l/2)r where / and m are any integers. It can be 
shown that 




E 




(12) 



n=— oo 




(13) 



and 




2 



z + ik. 



z — k. 



(14) 
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where d is the lattice constant and c is chosen to normalize z*) properly. ip\<_(z, z*) has 
regularly distributed zeros. For k = 0, the zeros of the two 8 3 functions coincide and they 
become double zeros. For k ^ 0, they shift by an amount — ik x and k y , respectively. 

Figure |3] shows energy (measured in units of Vq/tt) versus the flux ratio for the square 
lattice. Here the labels n and k have been suppressed, i.e, E nk values have been projected 
on to the energy axis. The energy spectrum is found to be self-similar, i.e. the same pattern 
is repeated with a different size and in a slightly distorted shape. With increasing flux ratio, 
the spectrum is squeezed since the amplitude of the matrix elements decreases with distance 
in a Gaussian way. 

In conclusion, a new method for the evaluation of the energy spectrum of a particle in 
the lowest Landau level moving in an arbitrary potential, is developed. The problem of 
Dirac delta scattering centers is solved for demonstration. The periodic case is found to 
have a self-similar spectrum. The method is very general since it can be applied not only to 
any distribution of impurities but any potential, as far as the particle is confined into the 
lowest Landau level. For the random distribution of scatterers |16|| , the calculated density of 
states is in perfect agreement with the analytic result of Brezin et al. [17 1. For a continuous 
potential, the matrix elements are labeled by continuous indices and hence the proposed 
approach may not be as practical as it is in the Dirac delta type potential case. 

This work was partially supported by the Scientific and Technical Research Council of 
Turkey under Grant No. TBAG- AY/123. One of the authors (Z.G.) is grateful to O. Tekman 
for helpful conversations. 
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FIGURES 

FIG. 1. Energy bands over one-quarter of the magnetic Brillouin zone for (a) a = 4/3 and (b) 
a = 10/7. The energy eigenvalues are plotted in units of Vq/tt. 

FIG. 2. Contour plots of |^k(r)| for a = 2. Zeros are marked with *. (a) k x = k y = 0, (b) 
k x = \/tt/2 and k y = 0, (c) k x = and k y = \/n/2, (d) k x = k y = \/tt/2. 

FIG. 3. Energy eigenvalues versus magnetic flux (in units of flux quantum) per unit cell for 
square lattice. 
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